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Abstract. We present adaptive optics high angular resolution (~0.1") spectro¬ 
scopic observations in the 3 pm region of eight well known Herbig Ae/Be stars 
with circumstellar disks. We detect the aromatic emission feature at 3.3 /im for 
four out of six of our objects with flared disks (HD 169142, HD 97048, HD 100453, 
HD 100546), some weaker additional features at 3.4 and 3.46 /im and nanodia¬ 
mond features at 3.43 and 3.53 (im in two of our flared object (HD 100546 and 
HD 97048 respectively). We also detect hydrogen recombination line at 3.74 /jm 
in practically all objects. The emission in the polycyclic aromatic hydrocarbons 
(PAHs) feature at 3.3 /im, additional and nanodiamond features in the 3.4-3.5 /tm 
region is, for the first time, spatially resolved in all the sources where the features 
are detected. The full-width at half-maximum sizes that we derive are typical of 
emission arising in a circumstellar disk. On the other hand, the continuum emis¬ 
sion is unresolved, with the exception of HD 97048 where it is marginally resolved. 
We compare the observed spatial distribution of the 3.3 /tm PAH feature and the 
adjacent continuum to the predictions of a disk model that includes transiently 
heated small grains and PAHs in addition to large grains in thermal equilibrium 


Habart et al 


2004a b The model predicts that, as observed, the 3.3 p m PAH 


emission feature is significantly broader than that of the adjacent continuum and 
that about 50% of its integrated intensity comes from a radius R < 30 AU. We 
find that the predicted brightness profiles reproduce very well the observed ones. 
This proves beyond doubt that the energetic 3.3 /im PAH emission feature takes 
its origin in the inner disk regions. 


Key words, stars: circumstellar matter - stars: pre-main sequence - ISM: dust, 
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1. Introduction 


Aromatic infrared emission bands (AIBs) attributed to PAHs have been detected in a 
number of pre-main sequence star of intermediate mass. According to the more recent 
statistics based on ISO s pectra, al most 60% o f the Herbig Ae/Be stars (hereafter HAeBe) 
show PAHs in emission (Acke fc van den Ancker 


In 


Habart et al 


200 


1 - 


(2004aJ), we have investigated the presence and properties of PAHs 


in circumstellar disks of HAeBe stars by comparing the predictions of disk models to 
observations with ISO and ground-based telescopes. We show that PAHs are present at 
the surface of disks, and absorb a significant fraction of the stellar radiation, similarly 
to PAHs in the interstellar medium. Consequently, PAHs are an important source of 
opacity in circumstellar disks and are likely to play a dominant role in the thermal 
budget ( via the photoelectric effect) and chemistry of the gas (via H 2 formation on dust 


surfaces, 


Habart et al 


2004b ). Moreover, the PAH emission is predicted to be extended, 


more that than of large grains at similar wavelengths. This is easily understood when 
considering the different excitation mechanism of the small transiently heated particles 
and of the large grains at thermal equilibrium. Because of the extension of the emission, 
which can be more easily resolved than that of the adjacent continuum with state-of-the- 
art adaptive optics (AO) systems and interferometers, one can use the small particles 
emission features to probe disk properties, such as the inclination on the line of sight of 
the inner disk. 

Recently, using AO high angular resolution spectroscopic observations of the HAeBe 
star HD 97048, we were able for the first tim e to spatially resolve the emission in the 


nanodiamond features at 3.43 and 3.53 //in (H aba rt et al 


2004c). We find that the 


diamond emission originates in the inner region (R < 15 AU) of the disk and show that 
the disk is seen pole-on, confirming previous indirect evidence. 

Here we present high angular resolution (0.1 ,/ ) spectroscopic observations in the 3 //m 
region obtained with NAOS-CONICA at the VLT of a sample of HAeBe stars. These 
observations allow us to spatially resolve the inner disk emission in the PAH emission 
feature at 3.3 //in due to C-H bond stretching vibrations, as well as in some additional 
and nanodiamond features in the 3.4-3.5 //m region. The paper is organized as follows. 
Sect. [2] presents our sample. Sect. |3] describes the observations and the data reduction. 
The results are presented in SectQ] In Sect. El we compare the observations to the model 
predictions. Sect. El summarizes our results. 


2. The sample 

Our sample consists of 8 well known HAeBe stars, all with evidence of a circumstellar 
disk and for which there are many observations at different wavelengths, including ISO 
spectroscopy. In Tabled! we report (i) the astrophysical parameters of each star (distance, 
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spectral type, effective temperature, luminosity, and y, i.e., far-ultraviolet (FUV) flux at 
10 AU from the star), (ii) the observational characteristics of the disk (evidence for 
a flared disk) and (iii) the dust spectral characteristics (presence of AIBs and silicate 
bands). We also give t he strength o f the 3 .3 u m aromatic emission fe ature measured 


from the ISO spectra of 


Meeus et al 


(2001) and 


Van Kerckhoven (2002). 


In this sample, we can distinguish three different groups (Meeus et al 


20Q1): the 


first (group la) is made of stars with evidence of a flared disk and strong or moderate 
aromatic and silicate bands; the second (group lb) is also characterised by flared disks 
and the presence of AIBs, but the silicate bands are absent; finally, in the third group 
(group II) there is no evidence of a flared disk, there are weak or no AIBs, but the silicate 
emission features are strong. 


3. Observations and data reduction 

Long slits spectra in the L-band (3.2-3.76 /im) were taken with the adaptive optics 
system NAOS-CONICA (NACO) at the VLT, using the visible 14x14 Shack-Hartmann 
wave-front sensor working at nearly 500 Hz. The 28" long slit had a width of 0"086 
which roughly corresponds to the diffraction limit while the pixel scale was ~0 / /0547 
and the spectral resolution ~700. For all objects, we took one slit position centered on 
the star. For one of our object, HD 97048, we took 7 slit positions: one centered on the 
star an d the others at 0"043, 0'/086, 0"129 and -0"043, -0"086, -0"129 off axis (see also 


Habart et al 


2004c ). The long slit was aligned in the north-south direction, except for 


resolved in scattered light (Augereau et al 


was aligned with the major 

, 2001 

Gradv et al.. 

2001) 


angle of ~160° measured east of north. For each object, in order to correct from the 
atmospheric and instrumental background, we employed standard chop/nod technique 
with a throw of ^9" in the north-south direction. The integration time per chop- or 
nod-positions was 1 min and the total integration time per source was about 30 min. 

A summary of the observations is given in Tabled The log of the observational condi¬ 
tions reports Fried parameters ro ranging from 7 to 25 cm (at 500 nm), average outer-scale 
£ 0 from 13 to 23 m, and a resulting average coherent energy (after AO correction and at 
2.2 fi m) from 25 to 60% on axis. The latter roughly gives the actual Strehl ratio charac¬ 
terizing the quality of our data. According to simultaneous measurements, the seeing (at 
500 nm) was between 0.6 and 1.5" consistent with the NACO values. Our observations 
clearly benefit from good seeing conditions and AO correction, and the achieved angular 
resolution was close to the diffraction limit. 


Data reduction was performed using a standard procedure for IR spectroscopic ob¬ 
servations. To remove the telluric features we used the observations of spectroscopic 
standard stars (listed in Table |2j taken for all sources immediately before or after the 
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Habartetal 


2004c '). These references stars provid- 


science target (except HD 97048, see 
ing a good match to the objects in terms of atmospheric turbulence conditions (ro, Cq, 
seeing) and AO correction quality were also used to estimate the point-spread function 
(PSF). 


4. Results 


4.1. Spectra 

Fig. H shows the spectra obtained for the sources in our sample integrated over the slit 
for distances from the star <1". When available, the figure also shows for comparison the 
spectra obtained with the Short Wavelength S pectrometer (SWS) on board ISO with a 
beam of ~14x20" around 3 pm and published in 


Meeus et al. 


(2001) and 


Van Kerckhoven 


l 2002). Both NACO and ISO spectra are normalized to the continuum emission adjacent 
to the 3.3 pm PAH feature. The wavelength coverage of the NACO spectra allows us to 
observe the 3.3 pm PAH feature, some additional and nanodiamond features in the 3.4- 
3.5 pm region, and two hydrogen recombination lines at 3.297 and 3.741 ^m, as discussed 
in the following. 


4.1.1. The 3.3 pm PAH feature 

In our NACO spectra, we clearly detect the 3.3 pm PAH feature for four out of six of 
our objects of group I, whereas we do not for the two objects of group II. 

For the foor group I objects showing the 3.3 pm feature (HD 169142, HD 97048, 
HD 100453, HD 100546), the feature is detected by both NACO and ISO. Comparing the 
NACO and ISO spectra, we find that the ratio of the emission peak to the continuum is 
roughly similar, i.e., that the larger ISO beam does not add more emission in the feature. 
This confirms that the PAH emission is confined to the disk. The exception is the ISO 
spectrum for HD 97048 which shows stronger PAH emission at 3.3 pm than the NACO 
one. In this case, one c an think that the surrounding nebula contributes to the emission 

(e-g-, 


van Boekel et al 


2004 ). 


The two group I objects where the 3.3 ^m feature is not detected are AB Aur and 
CQ Tau. In the case of AB Aur, the 3.3 pm feature is not detected by neither NACO 
nor ISO. On the ot her hand, PAH features at longer wavelengths (i . e., 6.2 , 7.7, 8.6 pm) 


have been detected (lYan Kerckhover , 


MB 


A eke Kr. van den Ancker 


200 41. One possible 


explanation of the non detection of the 3.3 pm feature is that in this object, which 
presents a particular strong near-IR continuum emission, the 3.3 pm feature disappears 
in contrast to the continuum, but other possibilities cannot be excluded, including a 
dominant contribution to the longer wavelength features from a residual envelope. 
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For CQ Tau, no ISO-SWS spectra around 3 fim are available. In this system, the non 
detection of the 3.3 fin i feature may be due to low excitation conditions due to the later 
spectral type (F2) of the central star. 

In the group II objects (HD 142666, HD 163296), an interesting case is HD 142666 
which shows in the ISO spectra weak 3.3 fin i feature and clear evidence of the 6.2 fim 
feature. In the NACO data, the signal-to-noise ratio is not sufficient to claim the presence 
or absence of a weak 3.3 fim feature. For HD 163296, there is certainly no strong 3.3 /rm 
feature but the data were noisy and we cannot say anything about a weak feature. 

In general, the NACO spectra indicate that PAH emission has size typical of 


previous studies based on ISO data (e.g., 

Meeus et al. 

Acke & van den Anckcr 

2004 

Habart et al. 

2004a 

)• 


2001 


Van Kerckhoven . 


2QQ3 


4.1.2. The features in the 3.4-3.5 fim region 


In two objects (HD 97048, HD 100546) the NACO spectra show, in addition to the 
aromatic band at 3.3 fim , some features in the 3.4-3.5 fin i region. 

HD 97048 shows the peculiar strong features peaking at 3. 43 and 3.53 u m, at¬ 
tributed to nanodiamond surface C-H s tretch i ng modes (e .g. 


Van Kerckhoven et a] 

Habart et al. 


200ii 


2m 


Sheu et al 


MB 


Jones et al. 



. As discussed in 


we see that the diamond components are in the NACO and 


ISO spectrum very similar in peak position, width and relative strength. Moreover, the 
peak/continuum ratio in the diamond bands is roughly similar. This suggests that nan¬ 
odiamond emission is confined to the disk. 

HD 100546 shows different, weaker features, peaking at 3.4 and 3.46 fim. In the NACO 
data, these features are clearly visible for distances from the star from ~ 0 . 2 " to 0 . 6 " (see 
Fig. E}. In t he IS O data, the shape of these bands is uncertain since the up and down 


scans differ (_Van Kerckhoven . 


200211 . These weak “satellite” bands have been observed 


ulae (e.g., 

Geballe et al. 

Van Kerckhoven 


2002 ). 


including reflection nebulae, planetary or p 

1985; 

Jourdain De Muizon et ah. 

1990a 

Q 


Joblin et al. 


1223 


200 21. They could be attr ibuted to aliphatic C-H stretches in methyl 


or ethyl side-groups attached to PAHs (e.g.. 


Joblin et al 


1996). 


4.1.3. The hydrogen emission lines 

Finally, the observed spectral interval contains two hydrogen emission lines, namely Pf 5 
close to the center of the broad PAH feature at 3.3 fim and Pf 7 at 3.74 fim. Pf 7 is 
visible in practically all objects and is stronger in the smaller aperture NACO data than 
in the ISO spectrum; its emission come mostly from the innermost regions of the system 
(see Sect. ITTal . 
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In the following, we use the observed Pf 7 intensity to estimate the contr ibut ion of Pf S 
to the 3.3 /im feature, using the theoretical intensity ratios given by Hummer & Storey 


198 71. For most of our objects, we find that the contribution of Pf S to the 3.3 /im 
feature is negligible (<4%). The exception is HD 100546 where Pf 7 is strong and Pf S 
seriously affects the observed shape of the 3.3 /im feature. For this object, we find that 
the contribution of Pf S to the 3.3 /im feature is strong close to the star and decreases 
very rapidly in the outer parts. For distances from the star <0.1", the spectrum at 3.3 /im 
is completely dominated by the contribution of Pf <5, we do no t see the PA H em issio n an d 
the measured Pf 7 /Pf S ratio is similar to that predicted by Hummer & Storev ( 1987 ). 
In the NACO spectrum for distances from the star < 1 ", the contribution of Pf S to the 
3.3 /im feature is about 40% ; this contribution is lower in the ISO data, about <7%, 
and the observed shape of the 3.3 /im feature is less affected. 


4.2. Brightness spatial distribution 

In this section, we present the brightness spatial distribution along the slit for the 3.3 pm 
aromatic feature, the features in the 3.4-3.5 pm region and the HI line at 3.74 pm.. We 
also show the brightness spatial distribution of the continuum emission adjacent to the 
3.3 pm PAH feature, which is due to warm large grains emission from the disk and/or 
stellar photospheric emission, as discussed in the following and in Sect. El 


4.2.1. Spatially extended 3.3 pm PAH feature 

Fig. El shows the intensity profile along the slit of the 3.3 pm PAH feature, the adjacent 
continuum and the observed reference star (the assumed PSF) for the 4 sources where 
we detect the 3.3 pm PAH feature. The 3.3 pm emission feature has been corrected for 
the contribution of Pf <5, as described in Sect. Era Errors on the measured full-width 
at half-maximum (FWHM) size values are dominated by systematic errors at the level 
of 20 %. 

In all objects the spatial extent of the PAH emission is significantly broader than that 
of the PSF. On the other hand, the continuum emissi on is unres olved, with the exception 
of HD 97048 where it is marginally resolved (see also 


Habart et al. 


2P_Q4d) . 


HD 169142. The PAH emission feature at 3.3 pm is extended with a FWHM of 0.3" 
(or 43 AU). The feature is clearly visible above the continuum and is strong in the inner 
regions (radius <0.3"), decreasing rapidly in the outer regions. For radius >0.3", the 
feature is weak and it is not possible to measure its brightness profile with sufficient 
accuracy. 


HD 97048. The 3.3 pm emission feature appears extended as in HD 169142. On one 
section of the slit and for distances from the star from 0 to 0.2", the NACO spectrum was 
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noisy near the 3.3 (im feature and it was not possible to measure its brightness profile with 
sufficient accuracy. As in HD 169142, the 3.3 /im feature is weak in the outer regions and 
given the uncertainties of the data, differences in the wings between the two objects are 
not significant. The continuum emission is spatially resolved in HD 97048 with a FHWM 
of 0.180" (or 32 AU); this indicates that in this system the continuum emission is mostly 
due to warm large grains emission from the disk. For HD 97048 for which we obtained 
NACO spectra on several slit positions, we were able to reconstruct two-dimensional 
images of the inner disk_in b oth the contin uum emission and the nanodiamond emission 
features (see Fig. 4, 


Habart et al 


2004(3) . The roughly circular shape of the emission 


maps confirms the hypothesis of a pole-on disk in HD 97049 put forward by 
(1986). 


The et al 


HD 100453. The 3.3 /mi feature emission is marginally resolved, with a FWHM of 0.18" 
(or 20 AU). In this object, the PAH emission feature decreases very rapidly in the outer 
zones. This could be due to low excitation conditions due to the later type (A9) of the 
central star. 

HD 100546. The PAH emission spatial distribution in HD 100546 differs in several 
points from that found in the other objects. Firstly, we detect the presence of a gap in 
the innermost region (R < 0.05-0.1" or 5-10 AU). In this region, as discussed in Sect. 
14.1.31 the spectrum at 3.3 /im is in fact completely dominated by the contribution of 
Pf S and we do not see the PAH emission. This gap is also found for the emission of 
the additional weaker features in the 3.4-3.5 /im region (see Sect. where there are 

no HI lines. It could be due to the lack of the carrier of the bands in the regions very 
near the star or to the fact that the emission of these features disappears completely in 
contrast to the continuum. But it could also result from the pres ence o f a “ ph ysic al hole” 


as suggested by the analysis of the ISO spectra of HD 100546 by 


Bouwman et al 


(2003). 


In order to explain the amount of mid-IR dust emission and the observed crystalline 


features, 


Bouwman et al 


2003) propose that the disk should have a hole and a wall 


at about 10 AU. In the NACO data there is no gap in the 3 /mi continuum emission, 
as expected in the case of a “physical hole”. However, we cannot rule out, on the basis 
of our data, that the stellar photosphere contributes a major fraction of the 3 /im flux. 
Higher spatial resolution data combined with an accurate analysis of the spectral energy 
distribution (SED) are necessary to shed light on this very interesting object. 

We also find that the 3.3 /tm emission feature is at larger distances stronger in 
HD 100546 than in the other objects. We can easily trace the 3.3 /im emission fea¬ 
ture to distances of ~ 0.5 arcsec or 50 AU from the star. Because of the presence of the 
gap, it is not possible to determine its FWHM size, but we note that our NACO data are 
consistent wit h a FWHM of ~0. 2" (20 AU) recently measured from ISAAC spectroscopic 


observations Geers et al. 


2QQJ. 
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4.2.2. Spatial extension of the features in the 3.4-3.5 /tm region 


We now look at the spatial emission distribution of the features in the 3.4-3.5 /mr region 
observed in two of our objects (HD 97048, HD 100546) and compare it to that of the 
3.3 /mz PAH feature. 

Fig- 01 shows the intensity profile of the strongest feature at 3.53 /jm attributed to 
diamonds, observed in HD 97048. The intensity profile of the other strong diamond 
feature at 3.43 /tm is similar. The derived FWHM is of 0.23" (41 AU) in both the 3.43 
and 3.53 /im features (see 


Habart et al 


2004c l. Also, by comparing the spatial extension 


of these features to that of the continuum and of the PAH at 3.3 /im, we find that the 
diamond emission is extended more than the continuum but less than that of the PAHs. 
This is roughly consistent with the idea that the diamond carrier of the 3 /mr features 
observed in HD 97048 have sizes intermediate between those of la rge grains a nd of PAHs 
as suggested by detailed analysis of the I SO spectrum (see e.g. 


2002 ; 


Sheuetal 


12002 : 


.Tones et al 


Van Kerckhoven et al 


200 41. Finally, we note that recent AO observations 


at the Subaru telescope of Elias 1, which shows the same diamond bands at 3.43 and 
3.53 /im as HD 97048, show similar results (M. Goto, private communications). The 
observed FWHM of the diamond features is similar to that derived here for HD 97048 
and the PAH emission is also found to be spatially more extended than that of the 
diamonds. 

Fig-EJshows the intensity profile of the strongest additional feature at 3.4 /im observed 
in HD 100546. As for the PAH 3.3 /im feature, there is evidence of a gap in the innermost 
regions. Also, the emission of the “satellite” feature is extended as the PAH emission. 
This is roughly consistent with the idea that these features are due to aliphatic C-H 
stretches in methyl or ethyl side-groups attached to PAHs. 


4.2.3. Spatial extension of the hydrogen emission lines 

Finally, we look at the spatial extension of the HI Pf 7 recombination line at 3.74 /tm. 
Fig. El shows the intensity profile of Pf 7 observed in two of our objects (HD 169142, 
HD 100546) where we clearly detected it. The HI emission line is not resolved. It originates 
from ionized material very close to the star, possibly the ionized gaseous disk or matter 
accreting from the inner disk onto the star. 


5. Comparison to model predictions 

In this section, we compare the observed spatial distribution of the 3.3 /mi PAH feature 
to the predictions of a disk mode l that inclu d es P AHs and very small grains in addition 
to large grains. These models iHaba/rt, et all LofMal reproduce the strength of the most 
commonly observed PAH features (i.e., 3.3, 6.2, 7.7, 11.3 /tm) in many HAeBe stars. 
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5.1. Description of the model 


In 


Habart et al. 


1 2004al) . we have computed models of disks heated by irradiation from 


the central star, in hydrostatic equilibrium in the vertical direction, with dust and gas 
well mixed (flared disks). We have used a 2-layer model to calculate the disk structure 
and implemented a ID radiative transfer code to compute the emerging spectrum. 

The models are appropriate for disks that are optically thick to the stellar radiation. 
This is generally the case for disks around pre-main-sequence stars, up to very large 
radii. For simplicity, we limit ourselves here to a disk with 0.1 M 0 , an inner radius 
of Ri n = 0.3 AU (close to the BG evaporation radius) and an outer radius of Rd=300 
AU. The surface density profile is E = Eo x ( r/Ri) _1 with Eo= 2 10 3 g cm -2 . 

The dust is a mixture of large grains in thermal equilibrium, transiently heated small 
grains and P AHs. T he big grains consist of graphite and silicates with optical constants 
from IDraine 1198511 : they have an MRN size distribution (n(a) oc a -3 ' 5 ) between a mini¬ 


mum and a maximum radius ai =0.01 pm. and a u =0.74 pm for silicates and cii =0.01 pm 
and a u =0.36 pm for graphite. To this grain population, we add a tail of the very small 
graphite particles with the same size distribution [rUa) oc a -3 ' 5 ) between a; =10 A and 
a u =100 A. Optical constants are from 


Draine 


For the population of PAHs, the absorption cross section is from 


Li & Draine (2001). 


In the ISM, PAHs are made up of a few tens up to a few hundreds of carbon atoms; for 
reasons of simplicity, we take in our model one PAH size, Nc = 100. The hydr ogen to 
carb on ratio is H/C = fn x ( 6/Nc) 0 ' 5 (case of compact symmetric PAHs, see 


Omont, 


198 61 with fn - the hydrogenation fraction of the molecule - equal to 1 (i.e., essentially 


fully hydrogenated PAHs). With respect to the charge, to keep the model simple we 
neglect that probably we have in disks a mixture of neutral and charged PAHs and assume 
that PAHs are neutral. Finally, we take into accont that in a strong FUV radiation field, 
photo-destruction of PAHs will occur during multiphoton events where PAHs absorb an 
energy of more than 21 eV in a time shorter than the cooling time. In our model, we 
assume that PAHs are destroyed when the probability that they are at a temperature 


T > T 


see 


(T P1 , nr , =2000 K) exceeds a critical value p eV ap (Pevap = 10 8 ) (for more details 


Habart et al 


2004alh The effect of these assumptions is briefly discussed in Sect. 15.41 
Finally, the silicate abundance in dust is [Si/H] = 3 x 10 -5 , and the total carbon 
abundance in dust is [C/H] = 2.22 x 10 -4 . Of this, 23% are in PAHs, 10% in VSGs and 
67% in large grains. 

The model predicts an infrared spectral energy distribution showing PAH features at 
3.3, 6.2, 7.7, and 11.3 pm clea rly visible above the cont inuum emission most ly due to 


warm large grains (see Fig. 3 in 


Habart et al. 


,2004a ). In 


Habart et al. 


(2004aj), we have 


compared the model results to observations from ISO and ground-based telescopes of 
some thirty HAeBe stars, including all our targets except CQ Tau. Taking into account 
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the uncertainty due to our disk description and assumed dust model, we concluded that 
the model can account for the strength of the most commonly observed PAH features 
(i.e., 3.3, 6.2, 7.7, 11.3 pm) in most of the objects with flared disks. 

For geometrically flat disks, PAH features are predicted to be very weak, even when 
PAHs with standard properties are present on the disk surface. This agrees with the 
obervations showing that in objects with no evidence of a flared disk (group II) no (or 
weak) PAH features have been detected. 


5.2. Predicted spatial distribution 


The model predicts that the PAH emission is much more spatially extended than the 
adjacent continuum emission, as shown in Fig. |H1 Fig. EJa) plots as a function of the 
radius the cunrmulative intensity of the PAH feature and continuum emission at 3.3 /im. 
The continuum reachs 50% of its intensity at a very small radius (about 2 AU) while 
the feature does so at larger radii (about 30 AU). This behaviour basically reflects the 
different excitation mechanism of PAHs, which are transiently heated by individual UV 
photons, and of the continuum, which is dominated by the emission of big grains, in 
thermal equilibrium with the radiation field. Only big grains located very near the star 
are hot enough to emit at 3 /im, whereas also PAHs farther away can be excited and emit 
in the 3.3 /im feature. However, also the 3.3 /im feature is expected to be stronger near 
the star since the emission of PAHs roughly scales with the intensity of the FUV radiation 
field. Moreover, among the PAH featu res, the 3.3 /trn one is the least extended spatially. 

( 2004a ). also for transiently heated particles the features 


As discussed in 


Habart et al 


at short wavelengths are in fact stronger in the inner part of the disk (where PAHs are 
hotter because multiphoton events occur where the radiation field is highest), decreasing 
rapidly in the outer cold part. On the other hand, the features at longer wavelengths are 
more extended. 

Fig-E^b) shows the surface brightness profiles of the 3.3 /tm feature and the adjacent 
continuum. As expected from Fig. Eta), the PAH emission feature is much more extended 
that the adjacent continuum. Finally, it is interesting to note that for the formalism 
adopted here, we find that very near the star (R <1-2 AU) PAHs are mostly photoe- 
vaporated, as shown by the inner gap seen in the 3.3 pm feature continuum subtraced 
emission profile (see Fig. 13b)). At R ~15 AU, about 50% of PAHs are photo-evaporated 
and for R > 20 AU this fraction goes to zero. 


5.3. Comparison with the observations 

Fig .[7|shows for each of our objects showing the 3.3 pm feature the intensity profile of the 
3.3 pm feature (continuum subtracted) obtained by convolving the computed intensity 
with the corresponding observed PSF. For each object, we have computed a model with 
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the suitable stellar parameters given in Tabled The disk and dust parameters are the 
same as in the previous discussed model shown in Fig. © and as described in Sect. EU 
The star/disk system is assumed to be seen face-on at the distance of the object. Fig. 
© also shows the predicted intensity profile of the adjacent continuum which is the sum 
of the computed disk continuum emission (mostly due to warm large grains) and of the 
photosphere stellar emission at 3 pm, each convolved with the observed PSF. The stellar 
contribution to the continuum emission at 3 pm is typically about 30% for our stars with 
spectral type of B9-A5, while it dominates for the later spectra type stars. The exact 
value of the stellar contribution depends not only on the stellar properties but also on 
disk parameters. Note, however, that the continuuum subtracted intensity profile of the 
PAH emission is independent of them. 

By comparing the predictions to the observations, we find that, as observed, the pre¬ 
dicted 3.3 pm PAH emission is significantly broader than that of the adjacent continuum 
and its predicted FWHM is very close to the observed one (differences < 20%, see Fig. 
0 . For some objects, there are intringuing differences between model predictions and 
observations in the wings, which are however generally within the uncertainties of the 
data. On the contrary, the adjacent continuum emission is not resolved or is only slighlty 
broader than the PSF. This also agrees with the observations. We take the general agree¬ 
ment between observations and predictions as a strong support for the physical pictures 
underlying our model. 

It must be emphasized that there are several complications which we have neglected. 
The most obvious is that we have assumed that PAHs can be characterized by a single 
size, hydrogenation and charge state. This is unlikely to be the case and is discussed in 
the next section. However, as discussed in the following, the predicted brightness profiles 
of the 3.3 pm feature for different PAHs properties are similar in the inner regions and 
show only significant differences at large radii, R > 40 AU. 


5.4. Discussion 


Our results allow us to locate the PAH emission in the disk, to show that the emission 
of the PAHs is spatially more extended than that of adjacent continuum and to prove 
beyond doubts that the energetic 3.3 pm PAH feature takes its origin in the inner disk 
regions. It is interesting to note that the spatial extension of the 3.3 pm feature is 
much smaller than those observed for PAH features at longer wavelength (i.e., 8.6, 11.3, 
12.7 pm), whi ch are exte nded on a scale of (several) 100 AU in HD 97048 and HD 100546 


van Boekel et al 


200 4). This is consistent with the model results, which showed how, 


among the PAH features in disks, the 3.3 pm one is the least extended spatially (see 
Sect. l5~2l) . 
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Our results have the interesting implication that PAHs survive photo-evaporation 
or coagulation processes even if these processes are in principle important in the inner 
disk, where the UV flux and the gas density are very high. Nevertheless, processes such 
as photo-evaporation must significantly affect the abundance of PAHs. For the formal¬ 
ism adopted here, we find that in the R < 20 AU region about 60% of PAHs (with 
Nc = 100) are photo-evaporated. But PAH evaporation is a complex process which 
increases rapidly with \ and for smaller PAHs and it is difficult to make quantitative 
estimates. Alternatively, it is possible that continuous replenishment of PAHs via subli¬ 
mation of icy mantles (in which interstellar PAHs may have condensed during the dense 
molecul ar c loud phase) or by accreting carbon atoms and/or ions from the gas (see, 


e.g. 


Allain et al 


199 61 maintain the PAHs population throughout the disk. This is a 


very interesting point and future detailed comparison with theoretical calculations will 
certainly add important information to the studies of the processes affecting the size of 
grains in disks. 

In the inner disk regions, processes such as ionization or de-liydrogenation are also 
expected to be important (because of the strong UV field) and to modify the PAHs 
properties. One expect PAH properties to vary as a function of the radius and depth in 
the disk. For example, PAHs are likely to be more positively ionized in the inner disk 
regions. Moreover, processes such as photo-evaporation or coagulation could affect the 
size of PAHs. Therefore, it is unlikely that PAHs are characterized by a single charge 
state, size or hydrogenation parameter, as assumed in our model. In order to get some 
insight into the specific PAH properties, one needs detailed com pariso n bet ween th e 


observations, not only of the intensity of the various features, as in 


Habart et al. 


i2004al) . 


but also of the brightness profiles at large radi i, R > 40 AU, where different PAHs have 
very different emission (e.g., Figs. 6 and 7 in 


Habart et al 


2004a b Unfortunatly, our 


NACO data in the wings are not of sufficient quality to be used in this context, but we 
want to emphasize that future higher sensitive spatially resolved spectroscopic data at 
~3 (im have the potential capability of constraining the PAH properties and of studying 
their evolution. 


6. Summary and conclusion 

We have presented NACO/VLT spectroscopic observations in the 3 /irn region of a sample 
of HAeBe stars. These data allow us, for the first time, to spatially resolve the inner disk 
emission in the bands of small transiently heated particles and thus provide a benchmark 
for testing the dust properties within the disk but also the inner disk properties. Our 
main results can be summarized as follows. 

1. We spatially resolve the 3.3 /irn PAH emission for all the sources where this feature 
is detected. We find that the PAH emission is significantly broader than that of the 
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adjacent continuum. On the other hand, the continuum emission is unresolved, with 
the exception of HD 97048 where it is marginally resolved. 

2. We spatially resolve the emission in the nanodiamond features at 3.43 and 3.53 //m in 
HD 97048 and find that the diamond emission is extended more than the continuum 
emission but less than the PAH emission. This is roughly consistent with the idea that 
the diamond carrier of the 3 pm features observed in HD 97048 have sizes intermediate 
between those of large grains and of PAHs, as suggested by detailed analysis o f the 


ISO spectrum (Van Kerckhoven et al 


2002 ; 


Sheu et al 


2002 : 


.Tones et al 


2004). We 


also detect weaker additional features at 3.4 and 3.46 pm in HD 100546, and find their 
emission is extended as the PAH emission. 

3. In HD 100546, we detect the presence of a gap in the innermost regions (R < 0.05- 
0.1" or 5-10 AU) in both the aromatic and weak additional features. This gap could 
be due to the lack of the carrier of the bands in the regions very near the star or to 
the fact that the emission of these features disappears completely in contrast to the 
continuum. But it could also result from the presence o f a “p hy sica l hol e” as s uggested 


by the analysis of the ISO spectra of HD 100546 by 


Bouwmanetal 


( 2003 ). Higher 


spatial resolution data combined with an accurate analysis of the SED are necessary 
to shed light on this very interesting object. 

4. We detect clearly the HI Pf 7 recombination line in practically all stars; the line is 
not spatially resolved, and probably originates from ionized material very close to the 
star, possibly the ionized gaseous disk or matter accreting from the inner disk onto 
the star. 

5. We compare the observed spatial distribution of the 3.3 pm PAH feature and the 
adjacent continuum to the predictions of a disk mod el that in clude s trans iently heated 


small grains and PAHs in addition to large grains fHabart et al 


i2004ah . The model 


predicts that, as observed, the 3.3 pm PAH emission feature is significantly broader 
than that of the adjacent continuum and that 50% of its integrated intensity comes 
from a radius R < 30 AU. We find a very good agreement of the predicted brightness 
profiles with our observations. This gives a strong support for the physical pictures 
underlying our model and show that the 3.3 pm PAH feature comes mostly from the 
inner disk regions. 

6 . Finally, our results have the interesting implication that PAHs survive photo¬ 
evaporation or coagulation processes even if these processes are in principle important 
in the inner disk, where the UV flux and the gas density are very high. Future detailed 
comparison with theoretical calculations will certainly add important information to 
the studies of the processes affecting the size of grains in disks. 
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Table 1 . Parameters of the targets. Evidence for a flared disk from spectral energy 
distribution modeling, infrared and millimeter interferometry observations. Presence of 
AIBs and silicate bands. Integrated strength of the aromatic 3.3 /im feature. 


Object 

d 

[Pc] 

Sp. Type 

T e ff 

[K] 

L* 

[L©] 

x“ 

Flared 

disk 

AIBs 

T h 

1 3.3 

[10~ 14 W/m 2 ] 

Sil. 

Group 

Ref. 

HD 169142 

145 

A5Ve 

10500 

32 

2.6 10 7 

V 

V 

1(0.2) 

- 

lb 

(1) 

HD 97048 

180 

B9-A0 

10000 

31 

2 10 7 

V 

V 

1.3(0.3) 

- 

lb 

(2) 

HD 100453 

114 

A9Ve 

7500 

9 

1.3 10 6 

V 

V 

1.3(0.2) 

- 

lb 

(1) 

HD 100546 

103 

B9Ve 

11000 

36 

3.6 10 7 

V 

V 

2.5(0.5) 

V 

la 

(1) 

AB Aur 

144 

B9/A0Ve 

9750 

47 

2.8 10 7 

V 

V 

< l c 

V 

la 

(1) 

CQ Tau 

100 

F2 

7500 

5 

7 10 s 

y / 

? 


? 

i 

(3) 

HD 142666 

116 

A8Ve 

8500 

11 

3 10 6 

no 

V 

0.3(0.2) 

V 

ii 

(1) 

HD 163296 

122 

A3Ve 

10500 

30 

2.4 10 7 

no 

- 


V 

ii 

(1) 


yj : detection, ?: possible detection, -: no detection. 

“ Far-ultraviolet (FUV, 6< hv <13.6 eV) flux at 10 AU from the s tar expressed in units of the 
average interstellar radiation field, 1.6 10~ 6 W m -2 lHahinJ . ll96iJ) . 

6 Integrated strenght (after continuum subtraction) and uncertainty (in between brackets) esti¬ 


mated from the ISO s pectr um of 


From 


Brooke et al. 


Meeus et al 


12001) and Van Kerckhoven ( 20021 1. 


1199311 


References: (1) 


erences therein; (3) 


Meeus et al. 


2001) and references therein; (2) IVan Kerckhovenl i2002fl and ref- 


Natta et al 


1 200111 . 
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Table 2. Summary of the observations. 


Object 

date 

ro 

[cm] 

Co 

M 

SR “ 

[%] 

seeing 

[arcsec] 

standard stars 

HD 169142 

2004/03/12 

7 

22 

25 

1.5 

HR 5494 

HD 97048 6 

2003/01/11 

10 

15 

40 

1.1 

HIP 53074, HIP 20440 

HD 100453 

2004/01/30 

14 

13 

55 

0.8 

HIP 062327 

HD 100546 

2004/01/31 

11 

16 

45 

1 

HIP 046283 

AB Aur 

2003/11/04 

11 

23 

30 

0.75 

GJ1055 

CQ Tau 

2003/12/02 

13 

20 

45 

0.6 

HIP025337 

HD 142666 

2004/03/07 

21 

22 

60 

0.75 

HIP 074752 

HD 163296 

2004/03/15 

25 

23 

60 

0.6 

HIP 085389 


a Resulting average coherent energy (after AO correction and at 2.2 /im) giving the actual 


Strehl ratio characterizing the quality of our data. 
b The data obtained for HD 97048 has been previously presented in 


Habart et al 


1 2004dl . 
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3.2 3.4 3.6 3.8 3.2 3.4 3.6 3.8 

A,(|im) A,(|im) 


Fig. 1. NACO spectra obtained for distances from the star <1" (thick solid lines). The 
thin lines show the ISO-SWS spectrum, obtained with a beam of 14x20". We also show 
as solid vertical lines the PAH feature position and as vertical dotted lines the position 
of Pf 7 line at 3.74 pm. Both NACO and ISO spectra are normalized to the continuum 
emission adjacent to the 3.3 pm PAH feature. For HD 163296, the data were noisy and 
it was difficult to remove the telluric features. 
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d (AU) 

-60 -30 0 30 60 


d (AU) 

-72 -36 0 36 72 




d (AU) 

-46 -23 0 23 46 


d (AU) 

-40 -20 0 20 40 




Fig. 2. Normalized spatial emission profiles of the 3.3 /zm PAH feature (continuum sub¬ 
tracted, solid lines) as function of the distance from the star. Dotted and dashed lines 
show the intensity profile of the adjacent continuum and of the PSF, respectively. The 
scale on the lower axis shows the distance from the star in arcsec, that on the upper 
axis the corresponding value in AU. For HD 97048, on one section of the slit and for 
distances from the star from 0 to 0.2", the NACO spectrum was noisy near the 3.3 /im 
feature and it was not possible to measure its brightness profile with sufficient accuracy. 
For HD 100546, we detect the presence of a gap in the innermost region (R < 0.05-0.1" 
or 5-10 AU); in this region, the spectrum at 3.3 /im is in fact completely dominated by 
the contribution of Pf 5 and we do not see the PAH emission (see text for more details). 
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d (AU) 

-72 -36 0 36 72 




Fig. 3. Left panel : Spectrum of HD 97048 in the 3.35-3.65 pm range obtained for dis¬ 
tances from the star <1". The dotted vertical lines show the strong diamond feature 
positions. Right panel: Normalized spatial emission profiles of the 3.3 pm PAH features 
(continuum subtracted, solid lines) and the 3.53 pm diamond feature (continuum sub¬ 
tracted, dotted-dashed line) as function of the distance from the star. The scale on the 
lower axis shows the distance from the star in arcsec, that on the upper axis the corre¬ 
sponding value in AU. 
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A, (pm) d (arcsec) 


Fig. 4. Left panel : Spectrum of HD 100546 in the 3.35-3.65 /jrn range obtained for 
distances from the star going from 0.2 to 0.6". The dotted vertical lines show the positions 
of the additional features. Right panel: Normalized spatial emission profiles of the 3.3 pm 
PAH features (continuum subtracted, solid lines) and the 3.4 /jm feature (continuum 
subtracted, dotted-dashed line) as function of the distance from the star. The scale on 
the lower axis shows the distance from the star in arcsec, that on the upper axis the 
corresponding value in AU. 



-0.4 -0.2 0 0.2 0.4 -0.4 -0.2 0 0.2 0.4 

d (arcsec) d (arcsec) 


Fig. 5. Normalized spatial emission profiles of the HI line at 3.74 pm (continuum sub¬ 
tracted, solid lines) as function of the distance from the star. Dashed lines show the 
intensity profile of the PSF. The scale on the lower axis shows the distance from the star 
in arcsec, that on the upper axis the corresponding value in AU. 
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(a) (b) 




Fig. 6. Disk model predictions. Panel (a): Cummulative intensity of the PAHs feature 
(continuum subtracted, solid line) and continuum (dotted line) at 3.3 /iin as function of 
the projected radius. Panel (b): Predicted surface brightness profiles at the peak of the 
3.3 /.mi feature (no continuum subtracted, dashed line) and in the adjacent continuum 
(dotted line). We also show as solid line the 3.3 /im feature continuum subtracted. The 
disk is heated by a typical HAe star with effective temperature T e g=10,500 K, luminosity 
L*=32 L 0 and mass M*=2.4 M 0 . The disk is flared, it has a mass of 0.1 M 0 , an inner 
radius of 0.3 AU and an outer radius of 300 AU. The star/disk system is seen face-on at 
a distance of 150 pc. 
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d (AU) 

-60 -30 0 30 60 


d (AU) 

-72 -36 0 36 72 
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Fig. 7. Comparison of model predictions and observations. Normalized spatial emission 
profiles of the 3.3 /im PAH feature (continuum subtracted, solid lines) and of the adjacent 
continuum (dotted lines) as function of the distance from the star. The thick and thin lines 
show the predicted and observed spatial emission profiles, respectively. For each object, 
the predicted emission profiles have been convolved with the corresponding observed 
PSF. The star/disk system is assumed to be seen face-on at the distance of the object. 
In each panel, the dashed line shows the normalized spatial emission profile of the PSF. 
The scale on the lower axis shows the distance from the star in arcsec, that on the upper 
axis the corresponding value in AU. 







